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When a classical black hole is perturbed, its relaxation is governed by a set of quasinormal modes 
with complex frequencies co — ujr + iuoi. We show that this behavior is the same as that of damped 
harmonic oscillators whose real frequencies are {uj\ + tjf )^", rather than simply ujr,. Since, for 
highly excited modes, wj ^ ojr, this observation changes drastically the physical understanding of 
the black hole spectrum, and forces a reexamination of various results in the literature. In particular, 
adapting a derivation by Hod, we find that the area of the horizon of a Schwarzschild black hole is 
quantized in units Ayl — 87r/|i, in contrast with the original result Aj4 = 41og(3)Zpi. 



Perturbations of black holes (BHs) vanish in time as a 
superposition of damped oscillations, of the form 



[asin(a;ijt) + b co's.{uj rI)] 



(1) 



with a spectrum of complex frequencies u = uir + iLOj. 
These quasinormal modes are of great importance in 
gravitational-wave astrophysics, and might be observed 
in existing or advanced gravitational- wave detectors. 
Furthermore, BHs are often used as a testing ground for 
ideas in quantum gravity, and their quasinormal modes 
are obvious candidates for an interpretation in terms of 
quantum levels. 

For Schwarzschild BHs, the quasinormal mode frequen- 
cies are labeled as w„i, where I is the angular momen- 
tum quantum number. For each I (with I > 2 for grav- 
itational perturbation), there is a countable infinity of 
quasinormal modes, labeled by the "overtone" number 
n, with n = 1, 2, . . .. In Fig. [l]wc show the frequencies of 
the I — 2 gravitational perturbations of a Schwarzschild 
BH of mass M: Imw„ grows monotonically with n, so 
the least damped mode corresponds to n = 1, and has 
2MImw ~ 0.1779 (we use units G = c=l). This is the 
mode that dominates the relaxation process. The next 
least damped mode is n — 2, with 2MImcL' ~ 0.5478, 
and so on. In contrast, the real part of co is not mono- 
tonic with n. It rather decreases at first, until it be- 
comes exactly zero for n = 9, and then starts growing 
again, reaching a constant asymptotic value. For large 
n the asymptotic behavior of the frequencies of gravi- 
tational perturbations is independent of I, and is given 

byaiiiaiaiiiii 



SnMuJn = \n3 + 2Tri{n+ -] +0 (n^^/^ 
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The pattern shown in Fig.[T]repeats for higher I. There 
is always a value fii of n such that, for n < hi, Re(w„^i) 
decreases with n, while above this critical value it raises 
again, up to the asymptotic value ln3/(87rA/) given by 
eq. 0. 

If we compare with the normal mode structure of fa- 
miliar classical systems, such as a vibrating rod, we have 
to admit that the structure displayed in Fig. [TJ and par- 
ticularly the "inverted branch" formed by the modes with 
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FIG. 1: Re{2MLu„i) against Im(2Mcj„i) for Z = 2 and n 
1,2,... 12, and for n = 20, 30, 40. Data taken from [J]. 



n < hi, is quite peculiar. In classical systems, the least 
damped mode is in general also the one with the lowest 
value of Re w, and typically Rew and Ima; both increase 
with n. In contrast, in Fig. [1] the least damped mode is 
the one with the highest possible value of Re lo and, for 
n < hi, Rew is a decreasing function of Imw! Even the 
"normal" branch n > hi is somewhat puzzling. Now Rew 
increases with n, which is more consistent with physical 
intuition, but still the fact that it saturates to a finite 
value is difficult to understand. In a normal macroscopic 
system, the underlying reason why, for large n, Ima;„ 
goes to infinity (and therefore these modes decay very 
fast), is that also Rea;„ diverges, so increasing n the 
wavelength Z„ — 27r/RecL;„ gets smaller and smaller, and 
finally becomes of the same order as the lattice spacing of 
the underlying atomic structure. At this point the per- 
turbation can no longer be sustained as a wave by the 
medium, and quickly disappears in the thermal agitation 
of the lattice nuclei. 

The quasinormal mode structure of Fig. [1] is no less 
puzzling if we attempt a semiclassical description and we 
intcrpretc it as the structure of excited levels of a quan- 
tum BH. In normal quantum systems, the levels with 
high excitation energy, E„ = hKeUn, are those that de- 
cay fast, first of all because, in a multipolc expansion, 
the decay width F grows with w (e.g. T ^ ui^ for a dipole 
transition and T ^ lu^ for a quadrupole transition) and. 



second, because they can decay into many different cfian- 
nels, i.e. into all the levels with lower excitation energy 
not forbidden by selection rules. So, again it is very sur- 
prising that, for n < fii, we have an inverted structure, 
where the lifetime of the state increases with its excita- 
tion energy. Quite puzzling is also the presence of a state 
with Rew„ = 0, and Inia;„ ^ (which exists for all /). 
So, the motivation of this work was to try to obtain a 
physical understanding of this level structure. 

To this purpose, we describe a quasinormal mode as a 
damped harmonic oscillator £,(t), governed by the equa- 
tion 



i + jd + ^U 



f{t)^ 



(3) 



where 70 is the damping constant, loq the proper fre- 
quency of the harmonic oscillator, and f(t) an external 
force per unit mass. Solving this equation in Fourier 
transform we get 
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div 



2tt [uj — uj+){lu — UJ-) 



(4) 



where uj± are the two roots of the equation w^ 



«7ow 



0, i.e. 



LU± — ± 



(7o/2)2 



.70 



(5) 



Consider the response to a Dirac delta perturbation, 
f{t) oc S{t), so f{cj) oc 1. For < < we can close the 
integration contour in eq. ^ in the lower half-plane and, 
since uj± both lie above the real axis, we get zero, as re- 
quired by causality. For i > we close the contour in the 
upper half-plane and we pick the residue of the two poles. 
So the result for £^{t) is a superposition of a term oscillat- 
ing as e'"+* and of a term oscillating as e"^"*. Therefore, 
the behavior ([T]) is reproduced by a damped harmonic 
oscillator, with the identifications 



Y = ^/ , v ^0 - (70/2)2 = ^^ 



which can be inverted to give 



ujo 



(6) 



(7) 



We see that the seemingly obvious identification loq = un 
only holds when 7o/2 <C wq, i.e. for very long-lived 
modes. For most BH quasinormal modes we are in the 
opposite limit; in particular, for highly excited modes, 
we have loi 3> (^_r, see Fig. [I] so wq — w/ rather than 
ujQ ~ LUjil If we model the BH perturbations in terms 
of a collection of damped harmonic degrees of freedom 
(which can be useful both classically, to have an intuitive 
physical picture of a BH as a whole, and in semiclas- 
sical quantum gravity, to get hints about the quantum 
structure of spacetime) the correct identification for the 
frequency of the equivalent harmonic oscillator is given 
by eq. Q, together with 70/2 = ujj. 
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FIG. 2: Re[2M{iuo)ni] for / = 2, against n. 



In terms of ljq the energy level structure of a BH be- 
comes physically very reasonable, and for ^ = 2 it is 
shown in Fig. [5] (a similar result holds for higher I). We 
see that the frequency ((jJo)n increases monotonically with 
the overtone number n. Recall that the damping coef- 
ficient 7o/2 is equal to loj, so also 7o/2 increases mono- 
tonically with n. Thus, in terms of the equivalent har- 
monic oscillators, the least damped mode, which is still 
the n = 1 mode, is also the one with the lowest value of 
luq, and the larger is {uja)n, the shorter is the lifetime, as 
we expected from physical intuition. 

For large n, using eq. ^ and introducing the Hawking 
temperature Th — h/{8TTM), eq. ([7]) can be written in a 
very suggestive form. 



h{uJo)7 



pI^ 



where 



mo = Tff In 3 , p„ = 2ttTh 



(8) 



(9) 



The expression for p„ is especially intriguing, since it cor- 
responds to a particle quantized with antiperiodic bound- 
ary conditions on a circle of length L = Ti/Th — SttA/ . 
It is also interesting to observe that the equal spacing of 
the levels for large n is just what would be expected from 
a description of the horizon in terms of an effective mem- 
brane [Ul . We can now reexamine some aspects of quan- 
tum BH physics, that have been previously discussed as- 
suming that the relevant frequencies were (a'/j)n, using 
(wo)^ instead. 

Area quantization. The idea that in quantum grav- 
ity the area of the BH horizon is quantized in units of 
Zpi (where lp\ is the Planck length) has a long history, 
that goes back to Bekenstein [12]. His result was that 
the area quantum of a Schwarzschild BH is AA = Stt^Ij. 
Hod [5|] found a similar quantization, but with a different 
numerical coefficient, using the properties of quasinormal 
modes of Schwarzschild BHs. Since for a Schwarzschild 
BH the horizon area A is related to the mass M by 
A = IQttNP, a change AM in the BH mass produces 



a change AA = 32ttMAM in the area. Hod considered 
a transition from an unexcited BH to a BH in a mode 
with n very large. He argued that for large n Bohr's cor- 
respondence principle should hold, so a semiclassical de- 
scription should be adequate even in the absence of a full 
theory of quantum gravity, and concluded from eq. ^ 
that the minimum quantum that can be absorbed in this 
transition has AM = hui = ?iln(3)/(87rAf). This gives 
AA = 41n(3)/p[ (recall that in units G = c = 1 we have 
Zpj = h). The numerical factor 41n(3) generated some ex- 
citement because of possible connections with loop quan- 
tum gravity [13| (see however \XM)- 

This conjecture is stimulating, but suffers from a num- 
ber of difficulties. First of all, further studies showed that 
the factor 4 In 3, that has its origin in ujji, see eq. ([2]), is 
not universal. For instance, we can consider a rotating 
BH with angular momentum per unit mass a. Comput- 
ing the asymptotic behavior of the quasinormal mode 
frequencies of gravitational perturbations, one finds that 
the large n limit and the limit a — » do not commute. If 
we first consider the asymptotic value for a Kerr BH and 
then we let a ^ 0, cjh does not reduce to In 3/(87rM ), but 
rather vanishes as a^^^ 0, [la, lla [l3 ■ This means that 
the area quantum becomes arbitrarily small if we give 
to the BH an infinitesimal rotation. Similarly, study- 
ing charged BHs, one finds that ujh changes discontinu- 
ously if we interchange the limits Q —^ and n — > oo. 
Furthermore, the study of generic spin-j perturbations 
revealed that the leading asymptotic value of the quasi- 
normal mode frequencies is given by [^ 



^SttMui 



(1 + 2cos7rj) . 



(10) 



For gravitational perturbations (j = 2) and for scalar 
perturbations (j = 0) the right-hand side of eq. (fTO|) is 
equal to —3, and we recover eq. ^. For vector per- 
turbation (j = 1), the right-hand side of eq. (flO)) is 
equal to +1, and we get SirMujn = 2mn, so the real 
part is now zero rather than In 3, and the correspond- 
ing quantum of area would also be zero. E quat ion (jlOp 
holds also for half- integer perturbations [la, ll9| ! iii this 
case the right-hand side of eq. (fTO|) is equal to —1, and 
SirMun = 2'Ki{n + 1/2), so again HeuJn vanishes asymp- 
totically. The conclusion is that the asymptotic value of 
Rewn (and also whether p„ in eq. ^ is proportional to 
n or to n -|- 1/2) depends on the spin of the perturbation, 
and is not an intrinsic property of the BH. A similar non- 
universal behavior was discussed in ^] in a large class 
of BH models, that in the (r», t) plane effectively reduce 
to 2-d dilaton gravity. In conclusion, the area quanti- 
zation determined by Hod's conjecture does not refiect 
any intrinsic property of the BH, and the would-be area 
quantum vanishes in various instances. 

Another criticism that can be raised to the above 
derivation is that one has considered only transitions 
from the ground state (i.e. a BH which is not excited) 
to a state with large n (or viceversa). However, it is also 
legitimate to consider transitions n —f n' where n and 
n' are both large. Bohr correspondence principle, that 



was advocated above, actually only holds for transitions 
where both n,n' ^ 1, so these are in fact the only tran- 
sitions that should be considered within the above logic. 
Now, if we use eq. © , we see that in a transition n ^ n' 
with n, n' ^ 1, Recj„ changes by 0(l/n^/^). This means 
that in these transitions the area changes by arbitrarily 
small amounts. So, even restricting to j = 2 perturba- 
tion of Schwarzschild BHs, the area quantization holds 
only for a transition from (or to) a BH in its fundamen- 
tal state, while transitions among excited levels do not 
obey it. 

All the above difficulties are removed when, in Hod's 
conjecture, we use (ajo)n rather than (ijjR)n- We consider 
a transition n -^ n — 1, with n large. Then (a;o)n — {'^i)n 
and eq. ^ gives an absorbed energy AM — ?i[(ajo)n — 
(wo)„-i] = ?i/(4M), so 



AA = 32ttMAM = 87r/|,i , 



(11) 



which coincides with the old Bekenstein result. At large n 
all other transitions require a larger energy; e.g. n —>■ n — 
2 takes away about twice the energy, since for large n the 
(wo)n are equally spaced. Even if we dare to extrapolate 
to low n, where semiclassical reasoning might go wrong, 
we still remain with a non- vanishing area quantum, of the 
order of STr/pj. As it is clear from Fig. [51 the transition 
from n = 2 to n = 1 is the one with the smallest possible 
jump. Using the values of uiji and ujj given in [1], we 
find {uJo)n=2 ~ (wo)n=i — 0.2/(4M), corresponding to 
A A ~ 0.2 (87r/pj), while the transition from n = 1 to an 
unexcited BH has AA ~ 1.5 (87rZ|,j). 

Contrary to what happens for ujr, the quantum of area 
obtained from the asymptotics of (wq)™ is an intrinsic 
property of Schwarzschild BHs: for large n the leading 
asymptotic behavior of wq is given by the 0{n) term in 
w/, and it does not depend on the spin content of the 
perturbation, as we see from eq. (jlO[) . Furthermore, in 
contrast to what happens to ojfl, for toi the limits a — > 
and n — > cx) commute, and similarly for the limits Q — > 
and n^ oo 0, [H, [3, [13 • The result ^ can therefore 
be consistently taken as an indication of a quantization 
of the area of the horizon of a Schwarzschild BH. (The 
generalization of these results to other spacetimes might 
however be non-trivial, see e.g. [2l|.) In this context, it is 
useful to remark that a gedanken experiment with black 
holes reveals the existence of a generalized uncertainty 
principle, which implies a minimum length of order Ipi 
,22j . and which fits very nicely with the above result. 

Entropy and microstates. If, for large n, the horizon 
area is quantized, with a quantum AA = alp^ (where for 
us a = 87r while for Hod a = 4 In 3), the total horizon 
area A must be of the form A — NAA — Nalp^, where 
N is an integer. Observe that N is not the same as the 
integer n that we used to label the quasinormal mode. 
Even for a BH in its ground state, n = 0, A^ is very large 
since it must account for the area of the unexcited BH, 
N = A/AA = 16TrM'^/{alli). The famous Bekenstein 
formula associates to the BH an entropy S — yl/(4/p[), 
and therefore at level N we expect that a BH should 



have a number of microstates g{N) given by g{N) ex 
exp{(AA)iV/(4;|,i)} = exp{aiV/4}. One might try to re- 
strict the possible values of a as follows [23, UM ■ One ad- 
mits the presence of a subleading term in the Bekenstein 
formula, S — A/(4Zp[) + const., and fixes the constant 
requiring that, for A^ = 1, there is only one microstate, 
g{N) = 1. This gives g{N) = exp{(a/4)(iV - 1)}. One 
then requires that g{N) be an integer. This restricts a 
to be of the form a = 41nfc, with k an integer. The 
value a = 4 In 3 is of this form, which is not the case for 
a = 8tt. 

However, a number of objections can be raised to this 
attempt to restrict a. First of all, in the semiclassical 
regime where our results could be trusted, N is very 
large, of order of A/lp-^, so g{N) is the exponential of a 
a very large number. Even if the number of microstates 
must be an integer, there is no hope that a semiclassical 
computation can reproduce this number with a precision 
of order one, necessary to distinguish an integer from a 
non-integer value. In fact, this does not happen even 
in classical textbook computations in statistical mechan- 
ics. Furthermore, the above expression for g{N) assumes 
that the same area quantum AA is valid from large N 
down to A'^ = 1, where our semiclassical approximation 
is certainly unjustified. Indeed, we see from eqs. ([8]) and 
© that the levels are equally spaced only in the limit of 
highly excited modes, otherwise there are deviations. 



Using our value of a in 5 
order in the large N limit. 



aN/4 we find, to leading 



5 = 27riV + 0(1) , 



(12) 



and g{N) ex exp{27rA^}. It is quite interesting to observe 
that eq. P^ agrees with the result found in Refs. 25, 26|, 
with apparently very different arguments. In these works, 
using the periodicity of the Euclidean BH solutions, it 
was found that the entropy is an adiabatic invariant, 
with a spectrum given, through Bohr's correspondence 
principle, precisely by eq. (|12p . This argument required 
only standard rules of quantum mechanics, but it was 
somewhat speculative in that the rules were applied in 
Euclidean time. 

On the other hand, the periodicity of the Euclidean so- 
lution also entered implicitly our arguments, since it is at 
the basis of the analytic computation of the asymptotic 
quasinormal modes frequencies, eq. PH]) . So it appears 
that the periodicity of the BH solutions in Euclidean 
time, beside providing a quick derivation of the value 
of Hawking temperature, is also at the origin of the area 
quantization law. 
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